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We calculate the electrically induced spin accumulation in diffusive systems due to both Rashba 
(with strength a) and Dresselhaus (with strength /3) spin-orbit interaction. Using a diffusion equa- 
tion approach we find that magnetoelectric effects disappear and that there is thus no spin accu- 
mulation when both interactions have the same strength, a = ±/3. In thermodynamically large 
systems, the finite spin accumulation predicted by Chaplik, Entin and Magarill, [Physica E 13, 744 
(2002)] and by Trushin and Schliemann [Phys. Rev. B 75, 155323 (2007)] is recovered an infinitesi- 
mally small distance away from the singular point a — ±/3. We show however that the singularity is 
broadened and that the suppression of spin accumulation becomes physically relevant (i) in finite- 
sized systems of size L, (ii) in the presence of a cubic Dresselhaus interaction of strength 7, or (iii) 
for finite frequency measurements. We obtain the parametric range over which the magnetoelectric 
effect is suppressed in these three instances as (i) \a\ — |/3| < 1/mL, (ii)|a| — |/3| ^ 7Pf, and (iii) 
l"^! ~ l/^l ~ \/ ^ /mppi with £ the elastic mean free path and pp the Fermi momentum. We attribute 
the absence of spin accumulation close to q = ±/? to the underlying U (1) symmetry. We illustrate 
and confirm our predictions numerically. 

PACS numbers: 72.25.Dc 85.75.-d, 75.80.-|-q 



I. INTRODUCTION 

Many recent works have explored mechanisms to gen- 
erate spin accumulations and spin currents by pass- 
ing electric currents through spin-orbit coupled elec- 
tronic systems. On the theoretical side, two related 
mechanisms have been proposed in disordered metals 
as alternatives to spin injection via ferromagnets or 
by optical means. They are current-induced trans- 
verse spin currents or voltages, a.k.a. the spin Hall ef- 
fect j^i^i^i^i^i^i^i^i^i^i^i^ and current-induced spin polar- 
ization (ClSP)^iiiii^ii^iilii^ The interplaybetween the 
two effects has been investigated in Ref. [I^- These ef- 
fects have been to some extent demonstrated cxperimen- 
tallyr^>^>^>^i^ and recent theoretical works have ex- 
tended them to include the mesoscopic regime, where 
fluctuations of both longitudinal and transverse spin cur- 
rents in mesoscopic ballistic and diffusive systems are be- 
ing investigatedi^i^i22ii^i2^ Most remarkably, it has been 
found that the same universality that applies to meso- 
scopic charge transport^ also applies to mesoscopic spin 
transport 1^ 

The main focus of these theoretical as well as experi- 
mental efforts is to use spin-orbit interactions (SOI) as a 
way to couple external electric fields to electronic spins, 
the hope being to generate, manipulate and/or mea- 



sure spin currents and accumulations by purely electrical 
means. SOI, however, also has the undesired effect of 
randomizing electronic spinsi^ This dichotomy theoreti- 
cally limits the use of SOI based magnetoelectric devices 
as components of information processors to the regime 
where the size of the device is much less than the spin 
relaxation length. A way to increase the spin relaxation 
length has been proposed in Ref. [13 for systems which 
exhibit SOI of both the Rashba^ 

= a{p^(Ty - Pydx) , (1) 

and the Dresselhaus type^ 

= PiPxO-x - PyO-y) , (2) 

where ax,y are Pauli matrices. When the two interactions 
have equal strength, a = ±(3, the SOI rotates electron 
spins around a single, fixed axis. The spin along this 
axis becomes conserved while spins along the perpen- 
dicular directions undergo a deterministic rotation that 
depends only on the starting and cndpoints of their tra- 
jectory. In particular, spins are not rotated along closed 
trajectories, therefore mesoscopic systems exhibit nega- 
tive magnetoresistance when a = ^^^ ^35,36,37 j weak 
localization and not weak antilocalization, just as if SOI 
were absent. An effective spin randomization still oc- 
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curs if the system is connected to many external trans- 
port channels, where uncertainties in the position of in- 
jection and exit translate into uncertainties in the spin 
rotation angle, unless injected electrons are prepared as 
spin-eigenstates of the Hamiltoniani^ In Ref. [33 (see, in 
particular, Eq. (7) therein) spatially periodic modes in 
diffusive systems have been first described for the case 
of equal strengths a = ±/3, with spatial period given 
by the spin orbit length. These modes arc long-lived 
for these particular SOI strengths (and in the absence 
of cubic SOI) and are thus referred to as persistent spin 
helix )'^^i^^i^° i.e. spin polarization waves with specific 
wavevectors {px,Py) = (4TOa,0). 

Charge currents flowing through spin-orbit coupled dif- 
fusive metals can generate finite spin accumulations.— 
This magnctoelectric effect achieves one of the main goals 
of spin-orbitronics - creating a steady-state, finite mag- 
netization solely by applying an external electric field. 
The direction of polarization depends on the direction 
of the electric field and on the spin-orbit interaction. An 
electric field in x-direction leads to an accumulation in y- 
or cc— direction for linear Rashba [Eq. ([I])] or Dressclhaus 
[Eq. (21)] interaction, respectively. The magnctoelectric 
effect in presence of both Rashba and Dresselhaus inter- 
action has been investigated in Refs.— which predicted 
that the CISP is given by the uncorrelated sum of the two 
accumulations generated by the Rashba and Dresselhaus 
SOI independently of one another. In particular, these 
predictions imply a finite accumulation at a = ±/3 7^ 0, 
whereas symmetry considerations (to be discussed be- 
low) require the vanishing of CISP at this point. This 
motivates us to revisit this issue. 

The purpose of this paper is twofold. First, we re- 
visit the Edelstein magnctoelectric effect in presence of 
both Rashba and Dresselhaus linear spin-orbit interac- 
tion. Contrarily to Refs i^^'^^ , we find that there is no 
CISP in any direction at a = ±/3. However, the spin ac- 
cumulation of Refs J^ii^ is recovered at an infinitesimally 
small distance away from the singular point a = ±/3 in 
infinite systems. Our second goal is therefore, and per- 
haps more importantly, to figure out to what extent phe- 
nomena occurring specifically at a = ±/3 are physically 
relevant. To that end, we consider three possible devia- 
tions from the treatment of magnctoelectric effects given 
in Refsii^iiiii^ in the form of (i) finite-size effects, (ii) the 
presence of a cubic Dresselhaus interaction 

H3D = liPyPlc^y - PxPIcTx) (3) 

which is always there whenever a linear Dresselhaus in- 
teraction is present, and (iii) an AC electric field. We 
find that spin accumulations arc suppressed over para- 
metric ranges given in each of these three instances by 

(i) |a| — |/3| ^ l/mL which depends only on the linear 
system size L, and not on the elastic mean free pat h i, 

(ii) |a| - |/3| < 7p|, and (iii) |a| - < hy/u/mpFe. 
There is a symmetry at a = ±/3 that is responsible 

for the vanishing of the magnctoelectric effect. In order 
to expose that symmetry, we first note that a linear SOI 



can be considered as a non-abelian SU(2) gauge field with 
components 

Ax ^ -m{acTy + (3a3:), Ay = m{aax + Pcjy), A^ ^ 0. 

(4) 

The corresponding field strength has only two nonzero 
components, 

Fxy = -Fyx = i[Ax,Ay] = -m^{a^ - p^)a, . (5) 

They vanish for a = ±(3. Alternatively, one can consider 
the rotated Hamitonian given below in Eq. (jlOp for a ~ P 
and perform the unitary transformation U = e*'^"='^/^ to 
obtain 

where H± = ^± 2apx + V. We thus see that the SU(2) 
gauge field reduces to two conventional U(I) gauge fields 
in the Hamiltonians H±. This U(l) field is a pure gauge 
field, implying vanishing spin conductance. To show this, 
one can for instance consider the linear response expres- 
sion for the spin conductance in a two-terminal meso- 
scopic sample^ 

G^= f dxdx'Tr[G^(x,x')J'G^(x',x)j;], (7) 

where the trace is over the spin degree of freedom, the 
integrals arc performed over cross-sections C^.j of the two 
leads connecting the system to external terminals and the 
current operators Ji = (iV; — A)/™, = {Jj,(Tf^}. We 

write G«'^(x,x') = c/«'^(x, x')e±*^ (''-'''), where 
is the retarded (advanced) Green's function of the system 
in the absence of SOI. For a = ±/3, one can gauge the SOI 
field out of the current operators via the transformation 

g«A.(x-x') j^g-»A.(x-x') ^ ^ .g^ 

TO 

which simultaneously gauges out the spin dependence of 
the Green's fimction in Eq. ([7]). We thus obtain (yu = 
x,y,z) 

G^^ f dxdx'V;5«(x, x')V,:.g^(x', x)Tr[a,J = .(9) 

■JC,,Cj 

It is straightforward to see that this gauge argument also 
applies to CISP, since the latter is given by a formula 
similar to Eq. ([7]), with the operator Jj replaced by a 
Pauli matrix. 

This article is organized as follows. In Section [III we 
use spin- and charge coupled diffusion equations to cal- 
culate the spin accumulation generated by a charge cur- 
rent flowing in a bulk diffusive sample with Rashba and 
Dresselhaus spin-orbit interactions. This approach al- 
lows us to consider spin polarization in a finite size sys- 
tem (Sec. Ill A[) . an AC external electric field (Sec. Ill Bp 
and in the presence of a cubic Dresselhaus interaction 
(Sec. Ill C[) . Section III presents numerical results on a 
tight-binding Hamiltonian confirming our analytical pre- 
dictions. A summary of our results and final comments 
arc given in the Conclusions section. 
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II. ELECTRICALLY INDUCED SPIN 
POLARIZATION NEAR a = ±13 

We consider a disordered 2DEG with non-interacting 
electrons of mass m and charge e. Choosing coordinates 
X, y and spin projections ax, cry along the crystal axes 
[110] and [110], respectively— the system is described by 
the Hamiltonian 

= i_ + n(p) • (T + ^/(X). (10) 

Zm 

Here, p = {px,Py,0) is the electron's momentum, cr = 
(ca;, CTy, CTz) is a vcctor of Pauli matrices (we later use 
cr° = 1), and 
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nip)= ^kjPjek ^ {-{a- P)py,ia + f3)px,0) (11) 

k.j=l 

is the internal field due to Rashba and linear Dresselhaus 
SOI given in Eqs. (P) and with strength a and /3, 
respectively. The disorder potential V is due to static, 
short-ranged and randomly distributed impurities lead- 
ing to a mean free path £ = ppr/m, where t is the elas- 
tic scattering time and pp the Fermi momentum. The 
interplay of disorder and SOI is characterized by dimcn- 
sionlcss parameters = 2apFT, = 2I3pft {h = 1) 
measuring the spin precession angle due to Rashba and 
Dresselhaus SOI between two consecutive scatterings at 
impurities. Our treatment presupposes ^a,/3 <SC 1, which 
ensures that spin distribution functions vary slowly ev- 
erywhere across the sample. 

The coupled spin and charge excitations of the 
Rashba/Dressclhaus spin-orbit coupled 2D electron gas 
obey the following integral equation (summation over 
doubly-occurring indices is assumed) 

/dr' 
^Tr[a^Gi(r, r')a,G^_^(r', r)]5.(r', lu), 

where >5'x,j,.z(r, w) and S'o(r, w) ~ n{r,uj) are the spin 
and charge distribution functions, respectively. We ob- 
tain diffusion equations in the presence of both Dressel- 
haus and Rashba SOI by gradient expanding this integral 
equation. For /3 = 0, these equations have been derived 
using diagrammatic perturbation theory^, kinetic equa- 
tions^ and quantum Boltzmann equation approaches^. 
For finite a and we obtain the same diffusion equa- 
tions as in Ref. [3^ which we rewrite here for convenience 

dtn = D\7^n + K^_^d,Sy~K^_^dySx, (12a) 

dtSx = DW^Sx-K^s-cdyn-KpxS.-TxSx, (12b) 

dtSy = DV^Sy + K^_^dxn~ KpyS, -TySy , (12c) 

dtS, = DV^S. + KydySy+K^dxSx-T.S,. (12d) 

Here the spin-charge couplings A"^!^^, precession cou- 
plings Kp'^ and spin relaxation rates Tx.y arc given by 
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For a homogeneous sample with a homogeneous charge 
current density, it is tempting to assume homogeneous 
spin accumulations and ignore all partial derivatives of 
5^ to obtain 

Sx oc -Kl^^TxdyTi = - (a - (3)Tdyn , (14a) 
Sy (X K^_^Tydxn = {a + f3)Tdxn , (14b) 
= 0. (14c) 

for the bulk Edelstein CISP. We note the cancellation 
of the potentially singular (a ± f3) factors in r^^y. For 
a — > ±/3 the spin-charge couplings go to zero but this be- 
havior seems to be cancelled by the diverging spin relax- 
ation time to give finite spin accumulations at a = ±/3. 
However the same approach for a set to /3 from the out- 
set produces vanishing spin accumulations. The main 
reason behind this inconsistency is that one spin relax- 
ation time of the system diverges as a ^ ±/3. However, 
in a real, finite-sized system, the spin relaxation time 
is bounded from above by the typical time to escape to 
the leads. This is so, because leads provide spin (and 
charge) relaxation, which for a — ±/3 becomes the dom- 
inant spin relaxation mechanism. Finite-sized effects are 
thus expected to induce a smooth crossover to zero CISP 
as a — > ±/3. In the next section, we show that this is 
indeed the case. 

A. Electrically induced spin polarization in finite 

systems 

We assume a rectangular sample with SOI attached 
to two external reservoirs defining the current direction, 
and bounded by vacuum otherwise. We obtain for the 
charge distribution function 

n{E) = 2(1 - x/L)T{E - eV) + {2x/L)T{E), (15) 

where J^{E) is the Fermi function. The appropriate 
boundary conditions are that the spin accumulations van- 
ish at the reservoirs and the normal component of the 
spin current vanishes at the hard wall boundarie a^^i^^ . 

Solving the diffusion equations we obtain the maxi- 
mum spin accumulation within the SOI region for an 
electric field along the a;-direction: 

Sy = WG(l-l/cosh(mi|a-/3|/?i)),(16a) 

(177 

^2DEG = (a + /3)r— . (16b) 
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For a field in the y-direction, one has the same behavior 
for instead of Sy, withla + /3| in the argument of the 
cosh and S'2deg = —{a — (3)Tdn/dy. Eq. (fTS]) shows that 
the Edelstein CISP goes smoothly to zero for a = ±/3, 
with the width of the crossover set solely by the system 
size, generating a singular behavior only as L — > oo. The 
size of the crossover region is in particular independent 
of the mean free path £, hence of the strength of the 
impurity potential, since in our regime, ^a,i3 ^ 1, the 
spin-orbit relaxation length is independent of disorder. 
Away from a = ±/3, one recovers the standard CISP 
S'2DEG predicted in RefsJIii^. The validity of Eq. dTT 
is illustrated numerically below in Fig. [2l 



B. AC-field induced spin polarization 

We next discuss the frequency dependence of CISP 
due to an AC electric field within the framework of the 
diff'usion equations (|12ap . For a = ±/3, this problem 
has already been addressed by Raiche\*^, and we revisit 
it briefly only for completeness. In an infinite system 
the polarization is spatially homogeneous such that all 
derivatives of S*^ in Eqs. (jl2b|) - (jl2dp vanish. The result- 
ing bulk polarization then satisfies 



{-iu + Tx)Sx 

{-iUJ + Ty)Sx 

{—iuj + r^)5'^ 



0. 



(17a) 
(17b) 
(17c) 



Further neglecting the influence of SOI on n one flnds 
from Eq. (fTSj) that Vn = —2veEj and thus 

2veEy{a - /3)Re[r^/(r^ - luj)] , (18a) 
-2veE^{a -h l3)Ke[Ty/{Ty - lu:)] , (18b) 
(18c) 



Sx — 

Sy = 

s, = 



for flnitc but small ujt <^ 1. This result has been found 
ini^di. As for flnitc-sizcd systems, we sec that both spin 
accumulations vanish at a = ±/3 and that the result of 
RefSi ^^i^^ is only recovered at \a\ — \P\ ^ -y/w/ mpFi. In 
the limit lot —>■ the polarization vanishes at the singular 
points only. 



C. CISP in presence of a cubic Dresselhaus 
interaction 

A linear Dresselhaus SOI, Eq. ([T|). is always accompa- 
nied by a cubic Dresselhaus interaction, Eq. whose 
strength might or might not be much weaker than that 
of the linear SOI. Because the presence of a cubic Dres- 
selhaus SOI breaks U(l) symmetry at a = ±f3, whose 
presence is crucial to the vanishing of the CISP, we in- 
vestigate in this paragraph the effect that a cubic Dres- 
selhaus SOI has on the CISP close to those points. 

If the cubic contributions are weak we still expect a 
suppression of the CISP at a = ±/3 and that the ad- 
ditional spin relaxation due to H^jj renders the point 
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FIG. 1: Upper panel: Spin polarization Sx^:^a)/S: 
E||[110] (dashed) and Sy,(„)/S'y,(„^o) (solid line) for E||[110] 
as a function of Rashba SOI a/P for = 2Ppft = 0.1 and 
^-y = 0.02. Lower panel: a-dependence of the normalized spin 
polarization Sy^(^a,-y)/Sy^(a,-y^o) for E||[110], = 2Ppft = 
0.1, and = 0.0, 0.01, 0.02, 0.03. 



a = ±/3 non-singular in the absence of boundary ef- 
fects and at zero frequency. In the coordinates chosen 
in Eq. (|10[) the cubic term in the SOI Hamiltonian is 



H3B = ^7 {pI - pI) iPx<yy - PyOx) ■ 



(19) 



incorporated into the diffusion 
The relevant relaxation rates F,, 



which has to be 

Eqs. (Il2b|-([i2d|. ....V„.X. ...C^C..^.. 

and spin-charge couplings K'^_^ have been calculated in 
Ref. |42 and Ref. |43, respectively. In our notation they 
are given by 



Fa-.n — 



2r 



4r 



^7 
16r 



(20a) 



4(a±/3)2 



2(/3±a)7P^ + -7'Pf 



j^x.y 

s — c 



(20b) 



where = 2jp^T and the upper (lower) sign applies to 
the X (y) component. In the presence of cubic SOI the 
relation K^l!^^ = rTy^xict T P), which led to the cancella- 
tion of divergent terms in Eq. (|14am4c() . no longer holds. 
The polarization is given by Eq. ([T7]), 



2,,ee,^,T-^K^_^E,, 



(21) 
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where now T^^y and K^^^ are given in Eq ((20b)) . and 
^zfii/ is the totally antisymmetric tensor of order three. 
The CISP is a rational function of £,a,p.-y Fig. [1] shows 
the behavior of Sx^y in the presence of weak cubic Dres- 
selhaus SOI {^.y = 2jp'pT <^ as a function of a/p. 
In this case, the polarization Sy does not vanish pre- 
cisely at a = (3 but shows a feature in the vicinity of 
this point. The minimum and maximum around the 
feature are at a = /?(! T ^7/(^/32^2)). The zeros are 
at a = /3(1 - and a = P{1 - ^y/^£.p). Thus 

we conclude that a weak cubic Dresselhaus interaction 
regularizes the singularity of the CISP around a = ±/?. 
The suppression of the CISP occurs over a width cx 7p| 
around a = ±/3. The predicted analytical dependences 
of on Rashba and Dresselhaus SOI strengths in Eqs. 
(Pna|) . (pOb)) and ^ may serve as guidance when at- 
tempting to tunc quantum wells to the symmetry points 
a = ztp and demonstrate the vanishing of the CISP due 
to linear SOI at this point. 



III. NUMERICAL SIMULATIONS 

We now perform quantum tranport simulations 
demonstrating the suppression of the CISP around the 
singular point a ~ ±/3 for finite size geometries. To 
this end we consider coherent electron transport in a dis- 
ordered quantum wire of width W with linear Rashba 
and Dresselhaus SOI. For the calculations we use a tight- 
binding version of the Hamiltonian (|10[) that we obtain 
from a discretization of the system on a square grid with 
lattice spacing a. The Hamiltonian is H = Hq + Hso with 

Ho = -ty^(4.gCg+a.tT + 4.gCg+i/.(T + h.c.) (22a) 

q.a 

Hso = ^[-(tR + iD)(cJ_^Cg+s,x-cJ_^Cq+£,|) (22b) 

+i(iR - tB)icl^Cq+y,l + cl^^Cq+y^f) + h.c] . 

Here ^ {cq^a) creates (annihilates) an electron with spin 
a =t or I in z-direction on site q = {qx,qy)- The vec- 
tors X and ij have length a and point in x and y di- 
rections, t ~ l/2ma^ denotes the hopping energy, while 
tR ^ a/2a and to = /3/2a are the Rashba and Dressel- 
haus SOI strength, respectively, in terms of which the 

spin-orbit lengths are given by £so = 7rat/tR_D- Wc fur- 
thermore include spin-independent disorder of Anderson 
type in the region of length L, where the on-site energies 
are randomly box-distributed with Uq G [—U/2,U/2]. 
The disorder strength U determines the elastic mean free 
path e « 48at^/'^y/E^/U^, which we tuned to values large 
enough that the system is not localized, but much smaller 
than the size of the disordered region in all our simula- 
tions 
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FIG. 2: (color online) Normalized spin accumulation 
S'j/Z&DEG as a function of a//3 for fixed /3/2a — to = 0.15t 
(giving ~ 21a), U = 2t (giving £ ~ 8.5a) and Fermi energy 
Ep = O.bt, for different linear system size W = L — 70a (red 
squares), 150a (blue diamonds) and 310a (grey circles). Data 
are averaged over 5000 disorder configurations. The solid lines 
are tlie tlieoretical prediction, Eq. (|16|l . with renormalized 
bulk spin accumulation and system size, S2DEG fet&DEG 
and L L^t with Sat ~ 0.84, L^t ~ 39.3a for L = 70a, 
Sat ^ 0.93, Lflt « 69.7a for L = 150a and Sm ^ 0.93, 
iflt ~ 117.1a for L — 310a. The electric current is in the 
direction x \\ [110]. 

We obtain the local electron and spin densities 

ji=-iTr[G<(9,<z)] , (23a) 
S^^ -iTv [a ^G<{q,q)] (23b) 

at site q by numerically computing the lesser Green func- 
tion G^{q,q). To this end we employ an efficient re- 
cursive lattice Green function method based on matrix- 
reordering algorithms as described in Rcf. |49|. Wc calcu- 
late averaged quantities {Si) and (n), over several thou- 
sands of disorder configurations and over a rectangular 
region in the center of the disordered part of the wire. 
We compare numerical data with the analytical predic- 
tion of Eq. ([TB|) . In Fig. [5] wc show the normalized, 
spatially averaged spin accumulation, (Sy) / S2I)EG- as a 
function of a//? varying the linear system size L. As ex- 
pected, we find complete suppression of {Sy) at a = /3, 
in agreement with Eq. (jl6p . Moreover, the pronounced 
dip around a = P becomes sharper and sharper as L 
increases, and the numerical data are in good qualita- 
tive agreement with the predicted line shape, Eq. (fTB|l . 
in particular, they have the same parametric dependence. 
The agreement becomes even quantitative if one normal- 
izes the system size and the bulk spin accumulation in 
Eq. Hn]), as is done in Fig. O Wc justify this normaliza- 
tion by the effective reduction of the spin-orbit interac- 
tion in confined systems with homogeneous SOI)^ and 
the fact that « 2.5£ is barely in the regime of va- 
lidity 4o > ^ of Eq. (Uni). This leads to smaller bulk 
spin accumulations and a longer spin relaxation length 
y than the case in which the conditions 
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FIG. 3: (color online) Disorder-averaged normalized spin ac- 
cumulation (Sy) /Sy ■,2DEG, with 5';,;2DEG ~ OT {du / dx) , as a 
function of a) the mean free path £ (for fixed width W = 50a) 
and b) the width W of the wire (for fixed U = 2t, £ 8.5a). 
The electric current is in the direction x \\ [100]. Differ- 
ent data sets corresponds to different values of /3/a = n/15, 
n = 10 (black circles), 11 (red), 12 (green), 13 (dark blue) and 
14 (light blue). In both panels, other parameters are fixed at 
tn = 0-15t, Ef = 0.5t, L = 40a and data have been averaged 
over 3000 disorder configurations. 



?a,/3 ^ 1 and Lg <^ L arc completely fulfilled, and quali- 
tatively explains the rcnormalization of the effective sys- 
tem length and the bulk spin accumulation. We also 
note that finite-sized effects lead to deviations from our 
estimates i sa iSat^/^y/Ep /U'^ for the elastic mean free 
path, and that numerical estimates based on the average 
inverse participation ratio^ of eigenstates systematically 
give a larger value for £ for which £/£^q ~ 0.55. 

According to Eq. pG]). the suppression of the CISP 
is independent of the strength of disorder / the elastic 
mean free path of the sample, as long as one stays in 
the diffusive regime. This prediction is supported by our 
numerical calculations. We find that the spin accumu- 
lation stays approximately constant with respect to the 
electronic mean free path. This is shown in Fig. [3^. In 
Fig. [SJd we moreover confirm that the CISP is indepen- 
dent of the width W of the rectangular SOI region for 
W > L. However, we expect a width dependence in the 
form of a reduction of the CISP upon reducing W , when 
D'yakonov-Pcrcl' spin relaxation^i begins to be reduced 
and finally suppressed due to the lateral confinement.— 



IV. CONCLUSIONS 

In this work we have studied the electrically in- 
duced and spin-orbit mediated spin accumulation in two- 
dimensional diffusive conductors with emphasis on finite- 
size and finite-frequency effects. In the thermodynamic 
limit of extended systems with (linear) Rashba and Dres- 
sclhaus SOI the Edelstein magnetoelectric effect gives rise 
to finite spin accumulation up to suppression at the sin- 
gular point |a| = \(3\. However, in many experimentally 
relevant systems, additional time, respectively, energy 
scales come into play, such as in tranport (i) through 
mesoscopic samples of finite size, (ii) in the AC regime 
and (iii) through samples with cubic Dresselhaus SOI. 
We have shown, both analytically and numerically, that 
in these situations the singularity in the spin accumula- 
tion at \a\ = \P\ is widened to a dip. This suppression 
of the spin accumulation over a finite Q;//3-range close to 
a ~ ±P may have interesting implications with regard 
to other phenonema based on the Dyakonov-Perel spin 
relaxation mechanism. As but one consequence, finite- 
size effects may render the spin-field-effect transistor pro- 
posed in Ref. [S^ for \a\ — \P\ effectively operative even 
if the two linear SOI arc not precisely equal. This is so, 
because the spin rotation along two different trajectories 
with the same endpoints remains the same, even away 
from |a| = |/3|, if the trajectories are not too long. This is 
rcfiected in the finite width \a\ - < 1/mL of the CISP 
lineshape given in Eq. ()16|) . Furthermore, given that spin 
helices also emerge from Eas. (jl2ap and ([T5)) '^^i^^ . we con- 
jecture that it is either finite-size effects or the presence of 
a cubic Dresselhaus SOI, or both, that render persistent 
spin helices excitable some distance away from a = ±/3, 
and thus experimentally observable. 

While the present analysis is based on diffusive charge 
carrier motion, it would be interesting to investigate bal- 
listic mesoscopic systems and see whether our results ap- 
ply there or if our analysis has to be extended. Work 
along these lines is in progress. 
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